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CTs . 1. Introduction 

m 

Dissipative hydrodynamics has been applied to relativistic heavy ion collisions to 
understand the viscous nature of the quark gluon plasma (QGP), in particular, to explain 
the experimental data of the collective flow parameter v 2 by including the effect of shear 
viscosity [T]. 

One of the objectives of these attempts is to determine 77/s, where 77 is the 
shear viscosity coefficient and s is the entropy density, by comparing hydrodynamical 
calculations of v 2 with the measured values. However, in addition to the uncertainties 
on the initial condition for hydrodynamics and the breakdown of the traditional 
Cooper- Frye freeze-out procedure in dissipative hydrodynamics (2j [3] , there are several 
difficulties in these attempts. One should note that these studies are mainly based on 
three assumptions; 1) the effects of other irreversible currents, such as the bulk viscous 
pressure and heat conductivity, are smaller than those of the shear stress tensor, 2) 
rj/s does not strongly depend on temperature and 3) the effects of shear viscosity on 
observables are mostly characterized by one transport coefficient, rj/s. 

In this work, we would like to discuss the validity of these assumptions. First, 77/s 
is in general a function of temperature. While the exact temperature dependence of 
this ratio is still unknown, its qualitative behavior is somehow well understood. As was 
shown in Ref. [5], 77/s exhibits a complex temperature dependence showing a minimum 
near the QCD phase transition. Thus, it is not clear that such a temperature dependence 
can be approximately replaced by a constant value. 
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Second, we remark that the effects of the shear stress tensor (or any other 
irreversible current) cannot be characterized only by rj. This might be the case for 
a fluid described by the relativistic Navier-Stokes (RNS) equation. However, as was 
shown in Refs. [5j |6], hydrodynamical equations which violate causality, such as the 
RNS equations, are unstable and cannot be used for any application. To solve the 
problem of acausality and instability, we have to introduce (at least) one more transport 
coefficient called the relaxation time tr, which represents the retardation effect in the 
formation of the irreversible current. Then it is natural to expect that the relaxation 
time should affect the behavior of Vi or any other observable. 

Because of the relaxation time, a large 7] does not necessarily induce a large shear 
stress tensor. To illustrate this, let us consider the retardation effect of an arbitrary 
irreversible current J which is induced by a thermodynamic force F. For the sake of 
simplicity, we consider a constant F. Then, the current is given by 

J = DF(l- e~ t/TR ), (1) 

where D and tr are the transport coefficients for this respective irreversible current. One 
can see that if tr is small enough compared to the typical scale of the hydrodynamic 
evolution, J is simply proportional to DF. However, if tr is not small, we have to 
consider the competition of D and tr to quantify the dissipative effects. 

As a matter of fact, the time scale of the hydrodynamic evolution in relativistic 
heavy-ion collisions is of the order of 10 fm. On the other hand, as is shown in Ref. [7], 
the relaxation time of the shear viscosity r n can also reach values around 10 fm. In this 
sense, we should not ignore the effect of r n in the analysis of v<i- 

In this paper, we calculate the effect of the shear and bulk viscosities on elliptic 
flow by taking a realistic parameterization of the shear and bulk viscous coefficients, rj 
and C) and the corresponding relaxation times, r n and rn. We argue that the behaviors 
close to ideal fluid observed at RHIC energies may be related to non-trivial temperature 
dependence of these transport coefficients. 



2. Model 



In this work, we will use the memory function method [TOl ITT] , which is one of 
the formulations of relativistic dissipative hydrodynamics. For a general metric, the 
conservation of energy and momentum is expressed as 

= 0, (2) 

where denotes the covariant derivative. We use the Landau definition of the 
four-fluid velocity u M in which case the energy-momentum tensor is decomposed as 
T^ v = eu^u u — (p + n) A MV + 7r^, where, e, p, U and 7r MV are the energy density, 
pressure, bulk viscous pressure and shear stress tensor, respectively. The projection 
operator is defined as usual, = g^ u — u^u v . 



The effect of shear and bulk viscosities on elliptic flow 



3 



In the memory function method, any irreversible current J is induced by a 
thermodynamic force F as 

t r u»D^J + (1 + t r D^)J = F (3) 

This equation is derived to take into account time retardation effects and the extensivity 
of hydrodynamic variables. One can regard this equation as the generalized version of 
the Maxwell-Cattaneo equation. It should be emphasized that the existence of the non 
linear term trD^u 1 ^ is essential to obtain a stable theory as is discussed in Ref. [TT] . 

By applying this equation to the cases of the shear stress tensor and the bulk viscous 
pressure, we obtain 

T U u a D a Tl + Tl = -(6-T n W, (5) 

where o~ pu = D <p u u> and 9 = D^u^ are the thermodynamic forces. We use the tradi- 
tional notation = A pua ^A al3 where = \ (A^ a A^ + A^A ra - §A^A a/3 ) is 
the symmetric traceless double projection operator. 

The above hydrodynamical equation are solved numerically in three spatial 
dimensions. We use the Smoothed Particle Hydrodynamics (SPH) method. See Refs. 
[T2"| ITTj for details. We use approximately 80000 SPH particles with a smoothing 
parameter h = 0.6 fm [12] . 

For simplicity, we use a factorized initial energy density profile into its longitudinal 
and transverse parts, as proposed by Ref. p3] . The transverse part is parameterized as 
in Ref. [2]. The longitudinal part is taken as in Ref. [13], with the choice of parameters 
VGaus = 0.8 and f]fi a t = 4. The initial value of the irreversible currents will be taken 
as zero and the impact parameter is always set as b = 7.5 fm. We use the equation of 
state (EoS) introduced in [2]. As for the freeze out, we used the ordinary Cooper- Frye 
formula [12]. However, we did not consider non-equilibrium corrections to the one- 
particle distribution function, because, as was pointed out in Refs. [2j [3] , this correction 
gives rise to unphysical behaviors and is not reliable. However, this will not affect the 
conclusions of this paper. 

As mentioned in the introduction, the temperature dependence of the transport 
coefficients is very important. The microscopic formula to calculate these quantities is 
given in Ref. [7]. However, the application of this formula is still limited to leading 
order perturbative calculations. Fortunately, in this approximation rj and ( become 
equivalent to the ones calculated in the Green-Kubo-Nakano formalism and we can use 
these results. In this paper, we use the lattice QCD calculations [8] for the QGP phase 
and we use the results from Ref. [9] for the hadron phase. These calculations for rj/s 
were fitted and the result is shown in Fig. [2]^a). We can see that r]/s shows a minimum 
around T ~ 180 MeV. For the bulk viscosity coefficient (/s, we use the parameterization 
showed in Ref. [2] with a critical temperature of T ~ 180 MeV. 

Similarly, in leading order calculations, the microscopic formula predicts a simple 
relation for the shear viscosity, r n = rj/P. For the bulk viscosity, there is still no result 
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Figure 1. Our parameterization of the shear viscosity rj/s(lett panel) and the shear 
relaxation time r n (right panel). 
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Figure 2. The elliptic flow V2(pt) calculated for different choices of 77/s (left panel) 
and calculated including the effects of bulk viscosity (right panel). 



and we simply assume r n = 9£/(e — 3P). Since t w (r n ) is linear with the transport 
coefficient 77 ((), the critical behavior of rj (() is transferee! to that of (rn). Thus, tv 
will have a minimum around the QCD phase transition while rn will display a maximum. 
In this sense, the complex temperature dependences of the relaxation times are induced 
by those of rj and (. 
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3. Results 

In Fig. Eta), we show ^(pr) for different choices of r]/s and neglecting bulk viscosity. 
These calculations were made for thermal pions with a freeze-out temperature of 
Tf = 130 MeV. The solid line corresponds to the case where rj/s has the temperature 
dependence shown in Fig. E^a) whereas the dashed, dashed-dotted and dotted lines 
correspond to constant values of rj/s of 0.08, 0.12 and 0.16, respectively. One can 
see that the elliptic flow calculated with the temperature dependent rj/s can be well 
reproduced by the one calculated with an effective constant shear viscosity coefficient 
Veff /s = 0.12. 

The quantity Tj e ff/ S is very close to the minimum of rj/s. This can be explained 
by the behavior of the relaxation time, shown in Fig. E^b). As mentioned in the 
introduction, r n characterizes the time scale for the velocity gradients to be converted to 
viscosity. Because the fluid expands with comparable time scale with r n , the magnitude 
of the shear stress tensor does not become large despite the large values of rj/s achieved 
throughout the hydrodynamical evolution. 

It should be noted that rj e ff/s is not bounded by the minimum of rj/s. Let us 
consider another set of parameters, r n = 3rj/P. In this case, the effective viscosity 
rj e ff/s becomes smaller than the minimum value of rj/s. This result shows that the 
value of rj e ff/s is not a direct measurement of rj/s. It can come from the combined 
effect of rj/s and r n , and the small values of rj e ff/s found at RHIC can just as well be 
a manifestation of the large relaxation time in the hadronic and QGP phases. Whether 
rj e ff/s can be used to quantify the value of viscosity at the phase transition or not, 
depends on the magnitude of the relaxation time in these regions. 

In Fig. E^b), we show ^(pt) considering the combined effects of bulk and shear 
viscosities, with realistic transport coefficients. The solid line corresponds to the 
calculation with only shear viscosity, the dashed line to the calculation with only 
bulk viscosity and the dashed-dotted line to the calculation with both bulk and shear 
viscosities. We see that the effect of bulk viscosity is rather small despite the large values 
of £/s near the phase transition. This can be explained as an effect of the relaxation 
time tq, again. Since r n also becomes large there, the values of (/s do not generate a 
large bulk viscous pressure. Since ( decreases very quickly outside the phase transition 
region, we cannot see considerable effects of the bulk viscosity. 

We remark that all the calculations showed in this work were done in three spatial 
dimensions. However, we were mainly interested in t> 2 at zero rapidity, where we 
confirmed that the Bjorken scaling ansatz is a good approximation. Thus, all the results 
showed here are approximately equivalent to (2+l)D hydrodynamical calculations with 
the Bjorken scaling for the longitudinal direction. This happened mainly because of our 
choice of initial condition, which displays a plauteu in energy density in the longitudinal 
direction. Studies with different initial conditions, specially in the longitudinal direction, 
will be done in future. 
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4. Concluding remarks 

We have shown that due to the interplay between the effects of viscosity and relaxation 

time, it is possible to fit the behavior of collective flow in terms of one unique effective 

shear viscosity r] e ff/s. However, its value depends critically on the evolution of the 

system. Specifically, the value is affected by the time span for which the system remains 

in different phases. For the LHC energies, for example, we would expect a different 

value for the effective viscosity since the system will stay longer in the quak gluon 

plasma phase. However, the details of such changes are hard to calculate due to the 

uncertainties on 77. We also remark that the uncertainties on the relaxation time are 

even larger since not so many effort have been applied to calculate the temperature 

dependence of this transport coefficient (unlike the rj/s case). 
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